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We study the density dependence of the decay ty° — > 77 using the Skyrme Lagrangian to describe 
simultaneously both the matter background and mesonic fluctuations. The classical ground state 
configuration has different chiral properties depending on the skyrmion density, which is reflected 
in the physical properties of pion fluctuating on top of the classical background. This leads to large 
suppression at high density of both photo-production from the neutral pion and the reverse process. 
The effective charges of are also discussed in the same framework. 
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I. INTRODUCTION 

At high temperature and/or density, the properties of 
hadrons are expected to change dramatically and under- 
standing these changes under extreme conditions is im- 
portant not only in nuclear and particle physics but also 
in many other related fields such as astrophysics. Data 
from high-energy heavy ion colliders, astronomical obser- 
vations on compact stars and some theoretical considera- 
tions suggest that the phase diagram of hadronic matter 
is far richer than the simple confinement/deconfinement 
picture seen in finite temperature lattice QCD simula- 
tions and include interesting QCD phases such as color 
superconductivity. Moreover effective theories can be de- 
rived for these extreme conditions, using macroscopic de- 
grees of freedom, by matching them to QCD at a scale 
close to the chiral scale A x ~ 4^/^ ~ 1 GeV. 

Chiral symmetry, which under normal conditions is 
spontaneously broken, is believed to be restored under 
such extreme conditions by virtue of its seeming to go 
hand-in- hand with confinement. The value of the quark 
condensate (qq) of QCD is an order parameter of this 
symmetry and is expected to decrease as the temperature 
and/or density of hadronic matter are increased. Since 
the pion is the Goldstone boson associated with sponta- 
neously broken chiral symmetry, the various patterns in 
which the symmetry is realised in QCD will be directly 
reflected in the in-medium properties of the pion, such as 
its mass m w and decay decay constant f n . 

These quantities have been the subject of previous 
studiespH2,liiB in the formalism adopted here. In those 
works, the Skyrme picture is used to describe both pions 
and dense baryonic matter. The basic strategy of the 
approach begins with the Skyrme conjecture that a soli- 
ton (skyrmion) of the meson Lagrangian can be taken 
as a baryon, so that dense baryonic matter can be ap- 
proximated as a system of infinitely many skyrmions. 



The pion is then incorporated as a fluctuation over this 
dense skyrmion matter. The chiral properties of the clas- 
sical background solution is directly reflected in the pion 
fluctuation properties, which may be interpreted as in- 
medium modifications. If we accept that the Skyrme 
model can be applied up to some density, its unique fea- 
ture of a unified meson-baryon description provides an 
interesting framework to investigate nonperturbatively 
the meson properties in dense baryonic matter. That 
is, we do not have to assume any density dependence 
of the in-medium parameters. We work with a single 
model Lagrangian whose parameters are fixed for mesons 
in free space. Only the classical ground state describing 
the dense skyrmion matter becomes highly density de- 
pendent and thus naturally in turn so do the fluctuating 
mesons on top of this dense background. There are how- 
ever a few drawbacks in the approach. Firstly, the low- 
est energy configurations for skyrmionic matter are only 
available thus far for a crystal structure we cannot 
yet describe the liquid structure of normal nuclear matter 
nor its behaviour at high temperature. Next, there are a 
few undetermined parameters in the model. They limit 
us to a qualitative understanding of the related physics. 
Leaving these drawbacks to further improvements, in this 
paper we apply the approach to the anomalous decay of 
the neutral pion into two photons, ir° — > 77 at finite 
density. Though this electromagnetic process may not 
be regarded as so important in the fireball phase of rela- 
tivistic heavy ion collisions, it is potentially relevant, for 
example, through 77 — > ir° — ► vv, to astrophysical phe- 
nomena such as core collapse supernovae and neutron 
stars 

The dependence of the process 7r° — > 77 on tempera 
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ture has been investigated by various authors 111 EH 
ITlj . It is well understood that though the anomaly, which 
drives the decay, is temperature independent, nonethe- 
less the amplitude does depend on T through the phase 
space or through the modification of the thermal quark 
propagators in evaluating explicitly the triangle diagram. 
In particular, [Tol 1 1 1| consider temperature effects in 
the WZW effective Lagrangian which incorporates the 
anomaly. The coefficient of the anomalous term corre- 
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sponding to 7r° — > 77 can pick up temperature depen- 
dence from loops 0, which is close in spirit to what will 
emerge in our study of density dependence below though 
for us the coefficient modifications come not from loops 
but at tree-level in the background field approach. 

Several works have studied also the density depen- 
dence of the process, with quite different results emerging 
IH 13 111 In [12 the decay is computed from the 
triangle diagram but keeping pion quantities such as mass 
and coupling fixed at their in vacuo values with the re- 
sult that the decay width increases with density. In a 
diagrammatic approach in the three-flavour NJL model 
is used with the result that the decay width decreases 
with density despite a corresponding increase in the pion 
mass. In |l4j the neutral pion decay width was computed 
by assuming its free space form from the anomaly and 
inputting two, quite different, model (such as Nambu- 
Jona-Lasinio) scenarios for the behaviour of the pion 
mass and decay constant at the phase transition point 
in temperature or density with correspondingly different 
results, increasing and decreasing, for the decay width 
dependence on density. 

In distinction to this, by using the Skyrme Lagrangian 
we obtain all in-medium quantities such as pion mass and 
decay constant from the same Lagrangian from which 
we compute the decay width. We stress that though 
our approach uses some of the machinery of the Skyrme 
approach, the most important aspect of the calculation is 
that the density dependence is obtained from the same 
effective chiral Lagrangian as that for free space using 
the fluctuation formalism. Other details of the Skyrme 
model do not play a significant role. 

In 0| , a similar process tt° — > 77 was analysed for the 
generalized pion fluctuations and the generalized photons 
in the color-flavour-locked phase using the correspond- 
ing Wess-Zumino-Witten term [0| for that phase. It 
was shown that the decay of the generalised pion is con- 
strained by geometry and vanishes at large density. The 
color-flavour-locked phase can be taken as the high den- 
sity limit of our approach, in the region where a hadronic 
description is no longer valid. 

The paper is structured as follows: the next section 
reviews the Skyrme Lagrangian and the properties of the 
dense skyrmion matter. In section III we consider meson 
fluctuations on this matter background and extract the 
in-medium modifications to the pion observables. There 
is a brief statement of conclusions in section IV. 



II. MODEL LAGRANGIAN AND DENSE 
SKYRMION MATTER 

We begin with a modified Skyrme model Lagrangian 
[l7| . which also incorporates the scale anomaly of QCD 
in terms of a scalar dilaton field: 

C = f -f (j-^j Tr(d^d»U)+£ sk 



+^ir(j-) Tr(Z7t + tf-2) + £ W zw (1) 

1 TO 2 f 2 

+^ X d»X ~ ((x//x) 4 (ln(x//x) - l) + l) 

where U = exp(zr-7r// 7r ) G SU(2) and x is the scalar dila- 
ton field. The parameters in Eq. correspond to phys- 
ical properties of the corresponding mesons: and / w 
(m x and f x ) are the pion (respectively, dilaton) mass and 
decay constant in free space. The Skyrme term £ s k is the 
higher derivative term introduced into the Lagrangian to 
stabilise the soliton solution of the Lagrangian, namely 

C sk = ^Tr([L„L,} 2 ), (2) 

where = (d^W (i? M = U\d^U)). Finally, the 
Wess-Zumino-Witten term £\yzw is necessary to break 
the symmetry of Eq.Q) under U — > U' which is not a 
genuine symmetry of QCD. The corresponding action can 
be written locally as [18j 

Swzw = - 2^2 / e^d 5 xTv (L„ ■ • • L a ) (3) 

in a five-dimensional space whose boundary is ordinary 
space and time. For U € SU(2), namely for two flavors, 
Swzw trivially vanishes (for three flavors this gives the 
hypothesised process KK — > mm). This will change 
when we couple to photons, as discussed below. 

The model Lagrangian has a few parameters which are 
fixed by the meson dynamics in baryon-free space. For 
the pions, we can fix the associated parameters to the 
empirical values as = 93 MeV and m, = 138 MeV. 
On the other hand, for the dilaton field, there are no 
available well-established empirical values for the mass 
and decay constant (or equivalently the vacuum expec- 
tation value). There have been a number of theoretical 
discussions on the field itself [l^. We use the values 
phenomenologically determined in nuclear matter stud- 
ies |20|. There is another unknown parameter which has 
not been determined by any directly associated exper- 
iment, the Skyrme parameter e, for which we use the 
conventionally used value |2ll |22| . Due to the ambigui- 
ties in these model parameters our study should be taken 
as having only qualitative relevance. 

The Lagrangian is invariant under global charge ro- 
tations, U — ► U + ie[Q,U], where e is a constant and 
Q = diag(|,— i) is the charge matrix for light quarks. 
The coupling of the pions and the photons can be incor- 
porated by promoting this to a local symmetry, 5U — > 
e(x)[Q, U}. This can be done by replacing the derivatives 
acting on the pion fields (not the neutral dilaton field) by 
covariant derivatives, <9 M — > = + ieA^ where the 
photon field A^(x) transforms as — > A^ — (l/e)d tJ _£ 
and e is the charge of the proton. For example, the cur- 
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rent algebra term for the pions is rewritten via 



Wess-Zumino-Witten term and describe the anomalous 
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(4) 



On the other hand, as is well-known, minimal substitu- 
tion does not work in gauging the Wess-Zumino action. 
The so-called trial and error Noether method 0] gives 
the gauged Wess-Zumino-Witten action 



wzw 



(U, = r wzw (C/) - e / d 4 xA»J 



d*xe» vaf} (d lt A u )A a 



"24tt 2 

xfr [Q 2 L (j + Q 2 Rp + QUQU^Lp] 



(5) 



where the anomalous part of the electromagnetic current 
of pions is defined as 



1 



■>,, - 487r 2 £ ^J d x ( 6 ) 
Tr [Q(L v L a Lp + R v R a R )] 

The normal part of the current coming from the current 
algebra term (J3J is 

■£ = *f (j) Tr MR^-L^] (7) 

and is conserved by the classical equations of motion of 
the fields. 

For meson dynamics in baryon free space, the vacuum 
solutions for the pions and scalar field are 



U vac = 1, 



A. v 



fx- 



(8) 



The fluctuations on top of this free-space vacuum are 
then incorporated via 

U = U V = exp(if • tt/U), and X = fx + X- (9) 

To lowest non-trivial order we obtain 

C = ^d^TTad^TTa - ±m 2 7r a 7r a 

+\d lx xd ll x - \rn\x 2 
+^eA"(a M ir+7r- - d^-ir+) 

+ ^°^ f ^ (io) 

12ir 2 f% 



Thus we generate the usual pion-photon interactions. 
Note that the normal electromagnetic current of pions 
leads to the isovector current associated with two charged 
pions. Correspondingly, the anomalous current gives the 
isoscalar current for the neutral and charged pions. The 
two terms with N c in Eq. (|10|l come from the gauged 
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"7T°. At tree level, it leads 



to the standard 7r u — * 77 decay width in vacuo 



N r e 2 



647T \l2n 2 f^ 



(11) 



On the other hand, the nonlinearity of the Lagrangian 
supports soliton solutions (skyrmions) carrying nontriv- 
ial topological winding numbers. With Skyrme's conjec- 
ture interpreting the winding number as baryon number 
we may simulate dense baryonic matter by using the me- 
son Lagrangian as a system made of many skyrmions. 
Let the lowest energy configuration for a given baryon 
number density be U (r) = no + if ■ n and Xo(f)- The 
Classical lowest energy state of the multi-skyrmion sys- 
tem is a crystal and there has been intensive work in 
the late 80's on a model containing only pions Of 
course, such a crystal structure is unrealistic for a multi- 
baryon system and a randomising of the multi-skyrmion 
arrangement is an open problem. There have been a 
few pioneering works attempting to incorporate statis- 
tical fluctuations (for example, see |23j), which can be 
adopted in our approach in principle. However, none of 
the properties we consider below will depend explicitly on 
the specific periodic properties of the crystal structure. 

Two well-separated skyrmions are in a lowest energy 
configuration when relatively rotated in isospin space 
about an axis perpendicular to the line joining their cen- 
tres. We thus consider as the lowest energy state of 
skyrmionic matter at relatively low density that for a face 
centered cubic (FCC) crystal where well localised single 
skyrmions are arranged on each lattice site such that the 
twelve nearest skyrmions are each oriented corresponding 
to the above lowest energy configuration for a skyrmion 
pair. At higher density, skyrmion tails start overlap- 
ping each other and the system undergoes a phase tran- 
sition to a more symmetric configuration, the so called 
the "half-skyrmion" cubic crystal. There, one half of the 
baryon number carried by the single skyrmion is con- 
centrated at an original FCC site where Uo = — 1 while 
the other is concentrated on the links where Uq = +1. 
Such a system has an additional symmetry with respect 
to Uq — > —Uq, which results in the vanishing of (Uo). 
This is often interpreted as a restoration of the chiral 
symmetry in the literature 0. More precisely though, it 
is only the average value of U over space that vanishes 
while the chiral circle still has a fixed radius f n . We call 
this a "pseudogap" phase to distinguish it from the gen- 
uine chiral symmetry restored phase for which the chiral 
circle shrinks to a point at U = 0. 

The dilaton field, when introduced in the Lagrangian 
Eq.Q to restore scale symmetry, effectively plays the 
role of a "radial" field for U and the restriction on the 
chiral radius is then relaxed. The pseudogap phase still 
remains as a transient process unless the dilaton mass is 
sufficiently small |3J. Shown in Fig.l is a typical numer- 
ical result for the average value of no and Xol fx m ^ ne 
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in Eq.Q. Expanding in U n we thus obtain the La- 
grangian for the in medium fields 

C = ^G ab {r)d^ a d^ b -^Ml{f)n 2 



+ 5^X^X-iM x 2 (f)x 2 



2 eC(r)A> i (d ll n + n- -8^-^+) 
N ° e D(r\-7r°r uua0 F F a 
ieN r 



(13) 



F(f)eV af, A„dvW + d a ir-df,n° 



FIG. 1: Average values of a and xl fx °f the lowest energy 
crystal configuration at a given baryon number density. 



lowest energy field configuration for a given baryon num- 
ber density. In the figure, the baryon number density 
is given in units of normal nuclear matter density po. 
It should be taken only as a rough guide since density 
dependence scales strongly with the model parameters, 
especially on the "Skyrme parameter" e and the scalar 
mass m x as p p oc e 3 and p c oc to 3 ,. One can see nev- 
ertheless that as the density increases the average value 
(no) drops quickly and almost reaches zero around the 
density p pi where the system enters the pseudogap phase 
in the half-skyrmion configuration. Compared to the chi- 
ral limit = where (no) vanishes exactly at p = p p 
(as illustrated in Fig.l by the dashed line), the transition 
happens rather smoothly and thus can be called an ap- 
proximate or quasi- pseudogap phase. As we increase the 
density further, the system remains in this approximate 
pseudogap phase for some range of p but the average 
value oix/fx continues slowly to decrease. At density p c , 
the {xl fx) 7^ phase and the (x/ fx) = phase have the 
same energy. Then, at density higher than p c , the latter 
comes to have lower energy and finally chiral symmetry is 
restored. To conclude this discussion, the dilaton in this 
approach is important for setting up a consistent mech- 
anism of chiral symmetry restoration within the Skyrme 
approach. In the mean field approach we will use in the 
following the dilaton will not play any further role beyond 
this enabling a correct framework for change of symmetry 
properties approaching the phase transition. 



where we have written only the terms corresponding to 
those in Ea. (|l(Jfl to emphasise the in-medium modifica- 
tions of the various contributions. Note that this does not 
alter the anomaly structure of the Lagrangian, which is 
fully respected by the Lagrangian Eq.JTJ. Only the coffi- 
cients of the terms corresponding to each process receive 
local effective corrections from the background matter 
through the potentials: 

G a b(r) = (xo/f x ) 2 i n l S ab + n a n b ), 
Mi(r)=ml(xo/f x ) 3 n , 

M 2 (r) = m x ( Xo //x) 2 [31n(xo// x ) + 1] 
+ (U/fx) 2 mn a ) 2 



-6m 2 (l-n )(xo// x )], 



C(f) = ( X v/fxY< 
D(r) = n\ 
F(r) = ri\ 



'?>■ 



In a mean field treatment of these potentials, we may 
take their spatial averages so that they reduce to con- 
stants, namely (G ab (r)) = GS ab , (M 2 x (r)) = M 2 X , 
(C(r),D(r),F(r)) = C,D,F. Then the Lagrangian can 
be rewritten simply as 



C = 



+ l -elA»{d^*+v*- -0 M 7r-7r*+) 



*tivF a fi 



(14) 



127T 2 /* 



^r-e^A.d^+d^-dpTT* 



III. FLUCTUATIONS ON TOP OF THE DENSE 
SKYRMION MATTER 

We now introduce mesonic fluctuations on the dense 
baryonic medium just described via 

U = x/U^Uox/lU, X = Xo + X (12) 



where we have carried out a wavefunction renormalisa- 
tion 7r* = \/ (G)ir a for the pion fields. All the physical 
parameters are evidently modified by the medium. Their 
density dependence is given by the following relations 

fllh = Vg, 

ml/m^ = x/Mj/G, 

e*Je = (C/G), 

(eN c )^ 3 /(eN c ) = F, 

(e 2 N c )* 2 /(e 2 N c ) = D. 
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FIG. 2: The pion vertices for the processes and the corre- 
sponding quark processes. 



We now see that three different kinds of effective elec- 
tric charge of pions appear depending on the different 
electromagnetic processes involved: that of charged pi- 
ons (denoted by e*), that appearing in the vertices for 
7T° — > 77 and thirdly that for the process 7 — > 7r + 7r°7r~. 
All three charges were simply a unit of electric charge 
in the Lagrangian Ea. (|f U|> for the pions in baryon free 
space. However, as illustrated in Fig. 2, all three charges 
originate from the electric charge of the quarks and thus 
how each term gets the charge factor e in its effective ver- 
tex is based on completely different detailed dynamics of 
the quarks. The baryonic matter influences each process 
in a different way so that each effective charge develops 
its own density dependence. 

Above all, the "effective" electic charge of the charged 
pions illustrated as Fig. 2 (a) comes from the effect of the 
surrounding medium just as in the case of electrons in 
condensed matter. As for the other processes, we also 
have to take into account not only the modification of the 
electric charge of the quarks but also the changes in the 
number of colors N* involved in the triangular or square 
anomaly diagram of Figs.3(b,c). The latter is expected to 
take place through the modifications of the quark propa- 
gator and the quark vacuum. As explicitly shown in |16| . 
for example, the evaluation of the anomalous triangle di- 
agram in the color-flavor-locked phase does not pick up 
the N c factor from the quark loop. 

Note that in Eqs.lO F/D = (ng)/(no + ra§) < 1 
so that {eN c )*^ 3 /(eN c ) JV 3 < (e 2 N c )*^ 2 /(eN c )^2, which 



implies that fewer colors would be effectively involved in 
the box diagram of process (c) than in the triangle dia- 
gram of process (b). 

Shown in Fig. 3 is the density dependence of the pa- 
rameters appearing in Eas. (|14|) . Quantities here are nor- 
malised with respect to their values in free space and all 
are reduced by the effect of the baryonic medium. Only 
the pion mass appears stable at low density correspond- 
ing to the symmetry broken phase. The other quantities 
show quite strong dependence on the medium density 




FIG. 3: Ratios of the in-medium quantities to the free ones. 
In order to illustrate the strong parameter dependence of p p , 
we present two results obtained with (a) e = 4.75 and (b) 
e = 6.0. The other parameters are fixed as f n =93 MeV, 
m^=138 MeV, / x =240 MeV and m x =720 MeV. 



even at low density. On the other hand, the pion mass 
drops quickly to zero in the pseudogap phase, while the 
other parameters exhibit plateau behaviour. In the chi- 
ral symmetry restored phase, all the parameters rapidly 
go to zero after p c . 

As a consequence of changes in the in-medium param- 
eters, the 7r c 
result 



77 decay width is modified, with the 



,*3 



7r°77 



\N c e 2 



64tt \ 12tt 2 /; 



(15) 



That this result is identical in structure to Eq. (|ll|l up to 
the replacements m w — * m* , /„■ — * /* and so on, is a 
consequence of the preservation of the chiral anomaly in 
the presence of dense matter. Thus for example the pion 
mass coming from the phase factor of the process is re- 
placed by the effective pion mass. In obtaining Ea. (ll5|l . 
we have used a very naive mean field approximation. The 
higher order contributions in the medium-generated po- 
tentials of Eas. l|14|l can be incorporated systematically. 
However, as shown in the previous works on the other 
properties of pion 0, 0| , after the resummation they lead 
only a minor corrections compared to the zeroth order 
values that are nothing but those obtained in the mean 
field approximation. 

Temperature effects could have been revealed, for ex- 
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ample, in [TlT ] from the pion loops, because the pion prop- 
agator at finite temperature differs from that at zero tem- 
perature by discrete Matsubara frequencies. One may 
extract a density dependence by evaluating similar pion 
loops by using the modified pion propagator with effec- 
tive pion mass. We expect, however, that the result is 
less important than the dominant nonperturbative one 
through direct in-medium modification of the effective 
parameters. 

Since the dilaton field is decoupled from the anoma- 
lous WZW term of the Lagrangian and is neutral, this 
particle cannot be involved directly in electromagnetic 
processes such as as 7r° — > 77. The presence of back- 
ground medium could generate a it— matter— x coupling, 
which can be absorbed into the resummation process like 
other terms higher order in the potentials Eas. (|14|) and 
leads to small corrections on top of the mean field result 
we have pursued here. However, as discussed in 0], the 
dilaton plays its most important role in the chiral sym- 
metry restoration of the background matter. Thus, the 
dilaton field is indirectly relevant to the process n° — ► 77, 
as mentioned above, by providing for a consistent context 
of chiral symmetry restoration within which the anoma- 
lous process can be studied. We shall return to this again 
below. 

Numerical results for Ea. (|f 51) are presented in Fig. 4. 
One can see that the neutral pion decay process is 
strongly supressed as the matter density increases. This 
comes from the reduction in the pion mass (solid line in 
the inset graph) and from the reduction in the strength 
of the corresponding vertex (dashed line), where the for- 
mer plays more a dominant role in the pseudogap phase 
and the latter in the symmetry broken phase. 

We now compare our results to other work specifically 
focussing on qualitative behaviour, whether the decay 
width is enhanced or suppressed. To this end, it is useful 
to state the three broad approaches within which tem- 
perature/density dependence is reflected in pion prop- 
erties and chiral symmetry restoration at the critical 
point. Phenomenologically these different approaches 
centre around different ways of realising that the ratio of 
m-nf-K in medium vs. in vacuo must be equal to the cor- 
responding ratio for m q (qq) (Gell-Mann-Oakes-Renner 
relation), which must in turn vanish if chiral symmetry 
is restored at the phase transition. Also relevant is that 
upon restoration of the symmetry, the pion and its chiral 
partner the a, here represented by the dilaton, must be- 
come degenerate. Thus we can label as approach A the 
case of the pion mass increasing to become degenerate 
with its chiral partner, while f v to approaches zero to 
fulfill the ratio of GMOR relations. Approach B would 
have both the pion and its chiral partner decreasing in 
mass to become eventually degenerate, precisely the case 
for our work with the pion and dilaton becoming degen- 
erate. No experimental evidence yet invalidates one of 
these against the other. In comparing our results with 
the other aproaches mentioned in the introduction the 
difference in the final result for the decay width will ul- 



timately depend on which approach (A or B) has been 
realised for chiral symmetry restoration. Furthermore, to 
compare with works that only include finite temperature 
effects we will follow the general consensus and take our 
density dependence as qualitatively indicative also of the 
behaviour as a function of temperature. 

Firstly, our decay width suppression disagrees with |T^ 
which however does not incorporate medium (tempera- 
ture or density) dependences in quantities such as the 
pion mass and the decay constant. To that extent [T^ 
corresponds neither to approach A or B above and in 
fact can only be relevant for low temperature or density 
far away from any transition to chiral symmetry restora- 
tion. Our result agrees with in the decay width, 
though in our calculation the pion mass decreases with 
density while |l3| have chiral symmetry restoration re- 
alised via approach A: the pion mass increases to that of 
the a meson to achieve symmetry restoration. Our ap- 
proach is closest to 0| in that the temperature/density 
dependence of the pion mass and the decay width are 
incorporated. In particular a vanishing pion mass at 
the phase transition point corresponds to case II of 0] 
and approach B delineated above. However, in [l4| the 
mass is almost constant in temperature and density until 
the phase transition and then either suddenly increasing 
(Case I) or decreasing (II). Consequently, the tempera- 
ture/density dependence of the pion decay constant plays 
the most important role, which makes the decay width 
enhanced over some range in both cases but near the 
phase transition dropping suddenly to zero in case II. 

In our approach the pion mass is almost stable in the 
chiral symmetry broken phase, then decreases steadily 
through the pseudogap phase to vanish at the transi- 
tion to the symmetry restored phase, thereby realising 
restoration via approach B, as above. On the other hand, 
both the effective pion decay constant /* and the ef- 
fective values for (N c e 2 )* decrease monotonically in the 
chiral symmetry broken phase but become stable in the 
pseudogap phase. If we do not take into account the 
change in (N c e 2 ), the decreasing pion decay constant 
alone would lead to an enhancement of the decay width 
(as in case II of 14] ) in the low density regime where the 
chiral symmetry is still broken. Since (N c e 2 ) 2 decreases 
faster than /*, the decay width becomes decreasing in 
that region. This in medium modification of (N c e 2 ) then 
is the key result distinguishing our result for the ir° — > 77 
decay width from others while realising chiral symme- 
try restoration via decreasing pion mass. Fig. 3 and 4 
illustrate the degree of sensitivity of our results on the 
Skyrme parameters. Either way, we see in the two curves 
of Fig. 4 that qualitatively the decay width suppression 
is robust against such variations. This result smoothly 
matches to the results of where the analogous pro- 
cess of generalised pion decay tt° — > 77 is also suppressed 
in the color-flavor-locking phase which would take place 
at higher density above p c . 
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FIG. 4: Ratios of the decay width of 7r° — » 77 to the free space 
value. The curves of the inset graph are the ratios (m* jm^) A 
and ((eN c /f w y/(eN c /f n )f. A gain the results obtained with 
two different values of the Skyrme parameter: e=4.75(solid 
lines) and e=6.0 (dashed lines) 



anomaly which drives this process. The key new result 
emerging from this Skyrme approach is that (N c e 2 )* can 
also be medium dependent in the hadronic phase just as 
it does in the CFL phase |l5j . 

We repeat that this study of ir° — > 77 decay at finite 
density using skyrmions is qualitative, due to the naivety 
in the crystal structure of the background baryons, the 
absence of a treatment of the Fermi statistics of the nu- 
clear matter and the absence of loop contributions. 

Our result would appear to be not so relevant to rela- 
tivistic heavy ion collisions, for example the observation 
of suppression of neutral pions created in central vs. pe- 
ripheral collisions as seen at PHENIX 24]. However, the 
scenario of low temperature and high density matter is 
particularly relevant to the context of compact stars. Ap- 
plications of these results to this problem are underway. 



IV. CONCLUSIONS 

The main result of this work is given in Fig. 4 which 
shows suppression of the process ir° — > 77 computed in 
a unified framework of the Skyrme model for pion fluc- 
tuations in a dense baryonic matter background. The 
in-medium dependences of the pion mass and decay con- 
stant arise in the same effective chiral theory as the 
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